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Gas Dynamics in the LINER Galaxy NGC 5005: 
Episodic Fueling of a Nuclear Disk 

Kazushi Sakainoto^'^'^, Andrew J. Baker^, and Nick Z. Scoville^ 

ABSTRACT 

0> . We report high-resolution CO (1-0) observations in the central 6 kpc (1') of the LINER 

Q^ ! galaxy NGC 5005 with the Owens Valley Radio Observatory millimeter array. Molecular 

gas is distributed in three components — a ring at a radius of about 3 kpc, a strong 
central condensation, and a stream to the northwest of the nucleus but inside the 3 kpc 
^ ■ ring. The ring shows systematic noncircular motions, with apparent inward velocities of 

~ 50 km s~^ on the minor axis. The central condensation is a disk of ~ 1 kpc radius with 
a central depression of ~ 50 pc radius. This disk has a molecular gas mass of ~ 2 x 10^ 
Mq\ it shows a steep velocity gradient and a velocity range (~ 750 km s~^) 30% larger 
than the velocity width of the rest of the galaxy. The stream between the 3 kpc ring and 
Csl ! the nuclear disk lies on a straight dust lane seen in the optical. If this material moves 

O I in the plane of the galaxy, it lies at a radius of ~ 1 kpc but has a velocity offset by 

up to ~ 150 km s~^ from galactic rotation. We suggest that an optically inconspicuous 
Q\ . stellar bar lying within the 3 kpc ring can explain the observed gas dynamics. This bar 

- ' is expected to connect the nuclear disk and the ring along the position angle of the 

northwest stream. A position-velocity cut in this direction reveals features which match 
the characteristic motions of gas in a barred potential. Our model indicates that gas in 
the northwest stream is on an Xi orbit at the bar's leading edge; it is falling into the 



o 

'^ \ nucleus with a large noncircular velocity, and will eventually contribute ~ 2 x 10^ Mq 

^ ' to the nuclear disk. If most of this material merges with the disk on its first passage 

of pericenter, the gas accretion rate during the collision will be ~ 50 M© yr~^. We 

^ . associate the disk with an inner 2:1 Lindblad resonance and attribute its large linewidth 

H \ to favorably oriented elliptical orbits rather than (necessarily) to a large central mass. 

The 3 kpc ring is likely an inner 4:1 Lindblad resonance ring — or a pair of tightly 

wound spiral arms — arising at the bar ends. Both scenarios can explain the apparent 

noncircular motions in the ring. The high rate of bar-driven inflow and the irregular 

appearance of the northwest stream suggest that a major fueling event is in progress in 

NGC 5005. Such episodic (rather than continuous) gas supply can regulate the triggering 

of starburst and accretion activity in galactic nuclei. 

Subject headings: galaxies: dynamics and kinematics — galaxies: ISM — galaxies: 
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Introduction 



Observational Data 



NGC 5005 is an Sbc galaxy hosting a low ioniza- 
tion nuclear emission-line region (LINER) at its cen- 
ter (Keel 1983). We list some of the system's global 
properties in Table |l|. Ho et al. (1997) propose that 
the nucleus is accretion-powered, on the basis of a 
broad (1650 km s^^ FWHM) component in a fit to 
the nuclear Ha line. Risaliti et al. (1999) suggest on 
the basis of X-ray observations that the active nu- 
cleus is fully obscured by an absorbing column of 
Nh > 10^"* cm^^. An alternative driver of nuclear 
emission may be ongoing OB star formation which 
follows a starburst ~ 10^ yr ago (Goodrich & Keel 
1986). On larger scales, NGC 5005 is classed as SAB 
(moderately barred) by de Vaucouleurs et al. (1991), 
although the location of the bar-like distortion is not 
apparent in optical images. The galaxy is thought 
to comprise a binary pair with NGG5033, a spiral of 
similar size and redshift at a projected separation of 
~ 250 kpc (for d = 21.3 Mpc: TuUy 1988). Neither 
galaxy displays distortions to its optical morphology, 
suggesting that there has been no significant recent 
interaction. 

— Table — 

We first observed the molecular gas in NGC 5005 in 
the course of a survey of 20 nearby spirals (Sakamoto 
et al. 1999a), using the Owens Valley Radio Obser- 
vatory (OVRO) millimeter array. We found a strong 
concentration of molecular gas (~ 10^ Mq) in the cen- 
tral kpc and an [apparently] incomplete ring of gas 
at a radius of 3 kpc. The nucleus revealed a very 
steep velocity gradient in CO — AT^ ~ 700 km s~^ 
within 10" or AV/Ar > 800 km s'^ kpc"^ — con- 
sistent with optical observations (Filippenko & Sar- 
gent 1985). The combination of the massive nuclear 
gas concentration, the large central velocity gradient, 
and the nearly circular ring motivated a more de- 
tailed study of this galaxy's gas dynamics and their 
implications for nuclear fueling. Our original observa- 
tions had relatively low sensitivity and resolution: the 
5" X 4" (0.5 kpc) beam only marginally resolved the 
^ 1 kpc structure in the nucleus. In this paper, we 
report high-resolution and high-sensitivity CO (1-0) 
observations of NGC 5005 using data from three con- 
figurations of the OVRO array. These observations 
include data from the new U configuration, which 
with baselines ^ 400 m permits ^ 1" resolution at 
this wavelength. 



2.1. OVRO CO(l 0) observations 

Our new aperture synthesis observations of the 
CO(l-O) emission in NGC 5005 were taken with the 
OVRO millimeter array between February and May 
1998, for a single pointing center at aigso = 13'^08™3?r7 
and ^1950 = -l-37°19'29'.'0. The array consists of six 
10.4 m telescopes equipped with SIS receivers. Digi- 
tal correlators were configured to cover 336 MHz (874 
km s~^) with 4 MHz (10.4 km s^^) resolution. Ana- 
log correlators recorded 1 GHz continuum data from 
the signal and image bands. For gain calibration, we 
observed the quasar 1308-1-326 every half-hour; 3C273 
was observed for passband calibration. The chopper- 
wheel method was used to correct for the atmospheric 
transmission. An absolute flux scale (uncertain to 
^ 15%) was established by comparing the apparent 
strengths of 1308-1-326 and Uranus. We adopt 1.5 Jy 
as the flux density for 1308+326 in the spring of 1998, 
an increase over the 1.0 Jy estimated for the origi- 
nal [survey] data from April 1997. Our observations 
were obtained in the L, H, and [new] U configurations; 
when combined with the preliminary data, we had the 
equivalent of five full transits with projected baselines 
from 10.4 m to 400 m. The visibility data were cali- 
brated with the OVRO MMA software (Scoville et al. 
1993) and mapped with the NRAO AIPS package. 

Three combinations of configuration and weighting 
were used to make maps (Table |^), to which we will 
refer according to their low (4"), medium (2'.'7), and 
high (1'.'5) resolutions. Figure n^ presents all three sets 
of zeroth and first moment maps, while Figure shows 
the low-resolution channel maps. The fraction of the 
total flux recovered by each of our maps is estimated 
from comparisons of our data with single-dish obser- 
vations from the FCRAO 14 m and the IRAM 30 m 
telescopes (Young et al. 1995; Braine et al. 1993). 
Each data cube was corrected for the primary beam 
response of the OVRO antennas (a Gaussian with 
HPBW — 60"), convolved with the single-dish beam 
(45" for FCRAO and 23" for IRAM), and sampled 
at the pointing centers of the single-dish observations 
to obtain values comparable to the single-dish fluxes. 
Our low and medium-resolution maps recover about 
74% of the total flux in the central 45" and 93% in the 
central 23". The high-resolution map recovers 50% of 
the total flux in 45" but 75% in the central 23", where 
there is a strong concentration of molecular gas. We 
made no correction to our data for missing flux. We 



did correct for primary beam attenuation before mea- 
suring flux from our maps, although this correction 
was not applied to the maps shown in the figures. 

— Table | — 

— Figures - || — 

2.2. Supplementary data 

For comparison purposes, we have obtained data 
for NGC 5005 at several other wavelengths from archives 
and the literature. Figure shows a Stromgren y 
band (A = 550 nm, AA = 20 nm) image obtained with 
the 2.5 m Isaac Newton Telescope at the Observato- 
rio del Roque de los Muchachos, Spain. Dust lanes 
in this image trace the spiral arms, which if trailing 
imply (since the east side of the galaxy is receding) 
that the northwest half of the galaxy is nearer to us. 

— Figure || — 

Figure ^ shows a WFPC2 image obtained with 
the Hubble Space Telescope (HST) .^ This image was 
taken in the F547M filter (A = 547 nm, A A = 49 nm). 
Cosmic ray hits were removed and images from the 
four chips of the camera were combined with IRAF 
and STSDAS. The absolute position of the HST im- 
age relies on the HST Guide Star Catalog and is accu- 
rate to about 0'.'5. The sharpness of the dust lane to 
the northwest of the nucleus supports the view that 
this side of the galaxy is nearer, and that the spiral 
arms are trailing. 

— Figure ^ — 

NGC 5005 has been observed with the VLA on 
several occasions. A compact central source is only 
marginally resolved at ~ 1" resolution. Table S 
lists the measurements of the position of the central 
source. We also present the mean position of the ra- 
dio source, which coincides well with the dynamical 



center determined from our CO data (§3.1) 
— Table | — 

3. Results 



Our CO maps show an oval ring of '^ V diame- 
ter and a strong central condensation, confirming our 
previous observations ( ^akamoto et al. 1999a ) . In ad- 
dition to these two structures, we find a "northwest 



stream" between the ring and the central peak with 
unusually large noncircular velocities. We describe 
each of these three components in order of increasing 
radius. 

3.1. The nuclear disk 

The medium and high-resolution maps of Figure 
|l| now clearly resolve the nuclear gas disk. The disk 
has an extent of 20" x 10" in the medium-resolution 
map, elongated along the isophotal major axis of the 
galaxy (p. a. — 65°) as determined on large scales. Its 
extent is smaller in the high-resolution map, proba- 
bly because of lower flux recovery and lower surface- 
brightness sensitivity. The emission feature ~ 6" 
northwest of the nucleus is kinematically distinct from 
the nuclear disk; we discuss it in the next section. 
The total CO flux of the nuclear disk is Sco = 310 
Jy km s~^, which corresponds to Mgas = 2.3 x 10^ for 
the Galactic conversion factor A^Ha/^CO = 3 x 10^° 



cm (K km s 



(Scoville et al. 1987; Solomon 



* Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space 
Telescope Science Institute. STScI is operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc. under 
NASA contract NAS 5-26555. 



et al. 1987; Strong et al. 1988 ) and a factor of 1.36 
for heavy elements. 

We fit the velocity field in the nucleus to a tilted 
ring model using the AIPS task GAL. The high- 
resolution data were used for the model input; the 
derived kinematic parameters are listed in Table ^ 
The inclination and major-axis position angle of the 
nuclear disk agree well with those from fits to the 
isophotcs of the galaxy as a whole (see Table |l|). The 
slight differences (less than 10° for both position an- 
gle and inclination) may be due to noncircular mo- 
tions in the nuclear disk, a warp of the galactic disk, 
an incorrect assumption for the intrinsic axial ratio 
in the isophote fit, or statistical errors in both mea- 
surements. We will adopt the large-scale isophotal 
values for inclination and position angle in the rest of 
our analysis. The inferred dynamical center coincides 
within 0"5 with the mean position of the nuclear ra- 
dio source in Table 0. The CO emission is largely ax- 
isymmetric about the dynamical center, with a slight 
(~ 1" or 100 pc) shift of the emission centroid to the 
northeast. 

— Table I — 

In the high-resolution map, CO (1-0) emission weak- 
ens at the dynamical center, due to absence of gas or 
higher excitation at the nucleus. The S-shaped con- 
tours in the isovelocity map near systemic velocity 
indicate noncircular motions in the nuclear disk. The 
noncircular component of the velocity is directed out- 



ward on the disk's minor axis. The full range of ve- 
locities in the nuclear disk (in the medium-resolution 
data) is 740 km s^^, or 830 km s~^ if corrected for 
the inclination of 62°. This is about 30% wider than 
the HI line width of the galaxy, W^n — 556 km s~^ 
(not corrected for inclination: |de Vaucouleurs 
19!! 1), as has been noticed previously (Braine 



et al. 

et al. 

19!! 3). The large linewidth is confined to the nuclear 
disk, as seen from a position-velocity cut along its 
major axis (Figure 0). The dynamical mass and the 
mass fraction of gas within radii of 5" (0.5 kpc) and 
30" (3.1 kpc) are given in Tabic g. The total mass in 
the central kpc (2 x 10^° Mq, derived assuming circu- 
lar rotation) almost certainly overestim ates the true 
dynamical mass, as we will discuss in §4T; but this 
estimate imposes an upper limit on any large mass 
concentration in the galactic center. Figure ^ shows 
that the nuclear CO emission coincides with a dark 
arc (on the near side of the galaxy) in the optical. 
The center of the nuclear disk falls on a conical op- 
tical feature whose wide opening angle (~ 120°) may 
be consistent with outflow out of the galactic plane. 

— Figure |^ — 

— Table I — 

— Figure g — 

3.2. The northvirest stream 

The emission from the northwest stream between 
the nuclear disk and the 3 kpc ring has anomalously 
high velocities which are most easily seen in channel 
maps. In Figure g, strong emission is seen 5" - 10" 
northwest of the galactic center in at least the velocity 
channels from 938 - 1062 km s~^ (if not more). The 
same structure is seen with better spatial and veloc- 
ity resolutions in the medium-resolution channel maps 
(Figure 0), most clearly from 979 - 1083 km s"^ The 
emission is anomalous in that it is redshifted with re- 
spect to the systemic velocity although it appears on 
the eastern, approaching side of the galaxy. The (line- 
of-sight) velocity offset between this stream and the 
galaxy rotation estimated from emission elsewhere is 
as large as ^ 150 km s~^. This feature is definitely 
real: it is seen in multiple subsets of our data at the 
same location and velocity. If it traces material mov- 
ing radially in the galactic plane, the implied motion 
is at 200 km s~^ (about half the free-fall velocity) di- 
rected inward. 

— Figure |^ — 

The size of the northwest stream in the high- 



resolution map is 5" x 3" (0.5 x 0.3 kpc), with its 
major axis running in an east-west direction. Due 
to missing flux, this map may underestimate the fea- 
ture's true extent — possibly closer to ~ 1 kpc as 
seen in the medium-resolution map. The CO flux 
(in the medium-resolution map) and the CO-derived 
mass of molecular gas are Sqo = 20 Jy km s~^ and 
-^^gas = 1.5 X 10^ Mq respectively. The centroid of 
the component is at '^ 6 " northwest (p. a. = —70°) 
of the dynamical center. The galactocentric radius of 
this centroid is 1.0 kpc if the stream is moving in the 
galactic plane. The velocity gradient in the stream 
also runs from east to west with higher velocities to 
the east; the total linewidth is roughly 100 km s~^. 
Figure Q reveals that the northwest stream coincides 
with a dust lane observed optically. This lane lies to 
the west of the nucleus and runs from the nuclear disk 
in the east to the 3 kpc ring in the west. The molec- 
ular gas of the northwest stream falls in the middle 
of the dust lane, where there is a dark feature in the 
optical image. 

3.3. The 3 kpc ring 

The low-resolution maps (Figures |l| and ||) show a 
CO ring at a radius of ^ 30" (3 kpc) as in our previous 
observations. In these maps, the 3 kpc ring is now 
found to be nearly continuous, completely encircling 
the center with few gaps. The ring is nearly circular 
when deprojected, as seen from comparing it to the 
reference ellipse (an appropriately projected circle) in 
the maps. Most of the ring emission lies within the 
annulus from 20" - 40" (2-4 kpc) , although the outer 
boundary is uncertain because of the primary beam 
attenuation. The total mass of molecular gas in this 
annulus is 2.2 x 10^ Mq. The mean and peak surface 
densities in this region are (Sgas) = 40 Mq pc~^ and 
^gas''^ = 200 Mq pc~^ respectively, after correction 
for inclination. 

The kinematics of the ring seen in the channel 
maps (Figure 0) are intriguing. Emission in each 
channel is not reflection-symmetric about the major 
axis of the galaxy (i.e. the major axis of the el- 
lipse). The asymmetry is such that the emission on 
the far/near (i.e. SE/NW) side of the galaxy is consis- 
tently shifted towards the receding/approaching (i.e. 
to the NE/SW) side of the minor axis. This is most 
obviously seen at systemic velocity in the 938 km s~^ 
channel. The emission centroids from the ring and 
the nucleus are indeed aligned, but the alignment dif- 
fers by 30° from the galaxy's minor axis. Attributing 



this offset to an error in the position angle determined 
from the galaxy's outer isophotes is unwise: it is hard 
to make an error as large as 30° given the large incli- 
nation and undisturbed morphology.^ The adopted 
isophotal position angle is also supported by the fact 
that the highest and lowest- velocity emission from the 
ring appears on the galaxy's major axis (see the 658, 
689, and 1218 km s"^ channels in Figure |). The 
kinematic major axis of the nuclear disk has p. a. « 
70°, also in close agreement with the isophotal ma- 
jor axis. We instead ascribe the asymmetries in the 
channel maps to systematic noncircular motions in 
the ring. On the minor axis, these motions appear 
inward, since the emission on the far/near side has 
lower/higher velocity than systemic. The consistency 
of the asymmetries from channel to channel suggests 
that the noncircular motion is systematic and exists 
across the ring, although it is not observable near the 
major axis. 

In optical images (Figures ^ and 0), the ring ap- 
pears as a dark arc on the near side of the galaxy. A 
dark spiral lane clearly emerges from this part of the 
ring and winds to the southwest at larger radius. A 
counterpart to this spiral feature on the far side of the 
galaxy would be harder to see, but there arc hints of 
optically bright patches emanating from the ring to 
the northeast. We do not detect CO emission from 
the spiral arms located outside the ring. 

4. Gas dynamics in NGC 5005 

The contrast in surface density between the nuclear 
molecular disk and the surrounding region within the 
3kpc ring is striking; it is also similar to the gas mor- 
phology often seen in the nuclei of barred galaxies 
(Sakamoto et al. 1999a). Beyond its appearance, the 
degree of central concentration in this galaxy is (nu- 
merically) one of the highest in the Sakamoto et al. 
(1999b) survey, whose results show that barred galax- 
ies in general have higher degrees of gas concentration 
than their unbarred counterparts. We propose that 
NGC 5005 contains a small bar within its 3 kpc ring, 
which builds up the central gas condensation and is 
itself traced by the northwest stream. We locate this 
bar in a roughly east-west direction so that the north- 
west stream can be produced by a shock along the 
dust lane at its leading edge. The central depression 
in the nuclear disk may result from concentration of 



gas towards X2 orbits at the inner 2:1 Lindblad res- 
onance of the bar, although it may also reflect out- 
flow or a variation in excitation. The 3 kpc ring likely 
arises at the end of the bar, due to gas accumulation 
in either an inner 4:1 Lindblad resonance region or a 
pair of tightly wound spiral arms. A schematic view 
of our interpretation is shown in Figure H. We present 
below our detailed cases for the existence of the bar^ 
and the origin of the 3 kpc ring, along with a brief 
discussion of the role played by these structures in 
fueling the nucleus. 

— Figure || — 

4.1. Bar-driven gas dynamics 

Figure shows a position-velocity cut at a posi- 
tion angle of 80°, roughly along the hypothesized bar. 
The cut is 5" wide and therefore includes the emis- 
sion from the northwest stream. The p-v diagram 
shows two distinct kinematic components. The com- 
ponent with the steeper velocity gradient arises in the 
nuclear disk; the component with the shallower gradi- 
ent originates in the northwest stream [at (-1-5", -1-100 
km s^^)] and parts of the 3 kpc ring. The overall 
"tilted X" pattern is characteristic of p-v diagrams 
in models of bar-driven gas dynamics: each kinematic 
component corresponds to a major family of noncircu- 
lar orbits supported by a barred potential. Such pat- 
terns (also described as "figure-of-eight") have been 
successfully used to infer the presence of bars in the 
Galaxy (Binney et al. 1991) and in edge-on external 
galaxies (e.g. Garcia-Burillo & Guclin 1995; Kuijken 
fc Merrifield 1995| ). Similar patterns have also been 
observed in a number of barred galaxies with moder- 
ate inchnations (e.g. M83: Handa et al. 1990; NGC 
7479: Laine et al. 1999; UGC 2855: Hiittemeister 
et al. 1999). 

— Figure || — 

Figure |0| presents a model of gas orbits in a barred 
potential which reproduces the main features in the 
position-velocity diagram for NGC 5005. We calcu- 
late periodic orbits using the damped-orbit model of 
Wada (1994) and the formulae of Sakamoto et al. 
(1999a: Appendix B) for the case in which two in- 
ner 2:1 Lindblad resonances (ILRs) exist. We adopt 



^ So far, no HI observations to determine the kinematic major 
axis have been reported for this galaxy. 



^Eskridge (1999) confirms the existence of a bar in near-infrared 
images of NGC 5005 obtained for the Ohio State University 
Galaxy Survey. In addition, K' snapshot images from the Palo- 
mar 60-inch telescope are consistent with a bar of length ~ 1' 
at a position angle of ^ 80° . 



a pattern speed for the bar which locates these reso- 
nances at radii of 1.1 and 2.1 (in dimensionless units 
for which corotation occurs at 10.0); this choice is mo- 
tivated by considerations which folfow in i |4.2| . The 
dominant orbits outside the outer ILR and inside the 
inner ILR are parallel to the bar and belong to the Xi 
family. Between the ILRs, orbits which are perpendic- 
ular to the bar and belong to the X2 family dominate. 
In the left panel of Figure fol the xi orbits are hor- 
izontally elongated. At the ILRs, the position angle 
of the orbits differs by 45° from the bar major axis 
because of a damping term in the model (Wada 1994; 
Sanders & Huntley 1976). We note that the stream- 
lines in a cloud-based model like this one can have 
sharp turns (giving rise to large velocity widths) at 
the leading edge of the bar (e.g. Lindblad & Lindblad 
1994). Such features have been attributed by other 
authors (e.g. Regan et al. 1999) entirely to the shocks 
which appear in full hydrodynamic simulations. Our 
ability to match observations with a simpler (cloud- 
based) approach gives us confidence that its use here 
is reasonable. 

— Figure |l^ — 

The right panel in Figure |l^ shows a position- 
velocity cut through our model along the bar major 
axis. Gas particles on X2 orbits are plotted as crosses, 
those on xi orbits at the leading edge of the bar are 
plotted as diamonds, and those on xi orbits at the 
trailing edge of the bar are plotted as dots. The Xi 
gas at the bar's leading edge and the separate X2 gas 
inside the outer ILR (i.e. the diamonds and crosses) 
form a "tilted X" pattern very similar to that seen in 
NGC 5005. Moreover, the gas on the model X2 orbits 
has a maximum rotational speed much faster than the 
circular rotation speed for this potential (shown as a 
solid line). This excess arises because the tangential 
velocity near the pericenter of an elliptical orbit is 
higher than the circular velocity. Because the major 
axis of the bar is close to the major axis of the galaxy 
(i.e. the line of nodes) in NGC 5005, material on X2 
orbits will reach pericenter while moving nearly along 
our line of sight. The bar model thus provides a nat- 
ural explanation for why the lincwidth of the nuclear 
disk exceeds that of the galaxy's full HI line; we do 
not need to appeal to a large central mass. 

Gas on xi orbits at the leading edge of the bar, 
producing the shallower component of the "tilted X" 
pattern in the model p-v diagram, appears to have a 
velocity less than circular for NGC 5005. This deficit 
arises (as above) because the bar and galaxy major 



axes are nearly (though not exactly) parallel: stream- 
ing along the bar will occur more in the plane of the 
sky than along our line of sight. In the context of 
this model, the northwest stream seen in NGC 5005 
is clearly identified as gas moving towards the nu- 
clear disk on an xi orbit at the leading edge of the 
bar. More tentatively, we suggest that the emission at 
(—20", -1-100 km s~^) in Figure may also trace the 
leading edge of the bar on an xi orbit. This feature 
(a possible "east stream") can be marginally distin- 
guished inside the 3 kpc ring in the medium-resolution 
map of Figure |l|. Little emission is seen in the regions 
of the p-v diagram corresponding to xi orbits at the 
trailing edge of the bar, but this is not entirely sur- 
prising. Many barred galaxies show strong concentra- 
tions of gas at the shocks which form in leading-edge 
dust lanes. The likely reason is that gas in a bar is 
directly channeled to a nuclear disk through xi or- 
bits at its leading edge. After colliding with the nu- 
clear disk near pericenter and entering onto X2 orbits, 
most of the gas is unavailable to flow past the nucleus 
and trace the trailing edge of the bar. This distinc- 
tion between the behaviors of the stellar and gaseous 
components is emphasized by Bureau & Athanassoula 
(1999) and Athanassoula & Bureau (1999) in their 
extensive studies of p-v diagrams for barred galax- 
ies. Our own model confirms that p-v diagrams can 
sometimes tell us not only whether a bar supports two 
families of orbits, but also whether the gas on the Xi 
orbits is at the bar's leading or trailing edge. 

4.2. The origin of the 3 kpc ring 

The 3 kpc ring has a nearly circular shape but 
clearly noncircular velocities (§||^). Figure ^ shows 
position-velocity cuts along the galaxy's [isophotal] 
major and minor axes. The minor-axis p-v diagram 
clearly shows that the 3 kpc ring has some degree of 
noncircular (i.e. radial) motion. The gas to the south 
(north) of the nucleus has negative (positive) veloci- 
ties with respect to the systemic velocity. The radial 
motion is directed inward with an amplitude of about 
50 km s~^ on the minor axis. We consider two possi- 
ble scenarios for the origin of the ring and its inward 
motions — that it constitutes a canonical resonance 
ring at the inner 4:1 Lindblad resonance, and alter- 
nately, that it is actually a pair of tightly wound spiral 
arms. Although we prefer the former hypothesis, the 
latter can not be excluded with our current data. 

An inner 4:1 Lindblad resonance region occurs 
where i7 — k/4 (calculated assuming a purely ax- 



isymmetric potential) is nearly equal to the pat- 
tern speed of the bar Qp (see reviews by Buta & 
Combes (1996) and Sellwood & Wilkinson (1993)); 
it is more widely known as the "ultraharmonic res- 
onance" (UHR), nomenclature which we will adopt 
here. An "inner" ring which encircles a bar and falls 
inside spiral arms is often identified with the UHR 
(e.g. Schwarz 1984a; Buta 1986). Can the UHR 
ring scenario explain the inward velocities we see in 
NGC 5005? The pattern speed chosen to produce Fig- 
ure ^places the model's UHR at a radius of 6.3; the 
ratio -RuhrZ-Roilr ~ 3 successfully matches the ra- 
tio between the radii of the 3 kpc ring and the nuclear 
disk. If the 3 kpc ring is truly a stable ring rather than 
a transient feature or a pair of spiral arms, its con- 
stituent clouds must follow closed orbits in the frame 
rotating with the bar. In this case, the velocities of 
material in the ring can be estimated from its shape 
(Buta 1986). Figure |l^ shows a deprojected version 
of the CO(l-O) emission from NGC 5005. The 3 kpc 
ring is roughly circular, but is slightly elongated par- 
allel to the bar with an axial ratio ^ 0.8, consistent 
with the average axial ratio of 0.82 observed for in- 
ner rings in barred galaxies (Schwarz 1984b). In a 
frame rotating with the bar, gas motion will be coun- 
terclockwise anywhere inside corotation, including at 
the putative UHR. Near the minor axis of the galaxy 
(i.e. on the vertical line in Figure O) , gas on an elon- 
gated orbit in the ring will therefore be on its way 
from apocenter to pericenter and will appear to have 
an inward velocity. 

— Figure |l^ — 

A more quantitative approach is taken in Figure 
03, in which line-of-sight velocities in the 3 kpc ring 
are plotted with respect to azimuthal angle. It is evi- 
dent that circular motion alone fails to reproduce the 
observations, leaving velocity residuals up to ~ 100 
km s""'^. In the UHR model, we can assume that non- 
circular velocities in an elongated ring will be sinu- 
soidal with a phase delay roughly equal to the posi- 
tion angle of the bar. The dashed line in Figure |lj 
shows an elliptical-orbit model generated for this as- 
sumption which produces reasonable agreement with 
the observations, including the 50 km s~^ noncircular 
velocity on the galaxy's minor axis (azimuth = ±90°). 
Although we do not attempt parameter fitting for this 
model, the UHR hypothesis is qualitatively consistent 
with the pattern of noncircular motions in the 3 kpc 
ring. In this scenario, the spiral arms outside the ring 
need not be tightly wound, and the arm which runs 



east (west) on large scales probably starts from the 
east (west) end of the bar. 

— Figure || — 

An alternate explanation for the 3 kpc ring is that 
it is not really a ring at all, but instead is a pair of 
tightly wound spiral arms. Streaming motions in a 
spiral density wave will be directed inward along the 
minor axis for all radii inside corotation. Numeri- 
cal simulations (e.g., Roberts 1969, and Roberts & 
Hausman 1984) as well as observations of a number 
of spiral galaxies (e.g. M81, Visser 1980; M51, Vo- 
gel et al. 1993) show inward motions on spiral arms 
near the minor axis inside corotation. The velocity- 
azimuth angle plot in Figure n3 supports one element 
of the tightly wound spiral arm picture. The residual 
velocities from circular motion are (though somewhat 
scattered) sinusoidal with a period of 27r rather than 
27r/3, as expected for any ringlike feature inside coro- 
tation which is formed by density-wave spiral arms 
(Canzian 1993). A cautionary note for this scenario 
is the sharply different pitch angle displayed by the 
spiral arms at larger galactocentric radii. We would 
at least expect the eastern (western) spiral arm to 
start from the west (east) end of the bar, if the spiral 
arms start from the bar's ends as is usually the case 
in barred galaxies. However, this picture may still 
necessitate a zero or (unrealistically) negative pitch 
angle for the "arms" comprising the ring (see Figure 
pi] ). Precise determination of the galaxy's inclination 
and position angle (e.g. with HI mapping) and large 
field-of-view color index maps would both help deter- 
mine whether the 3 kpc ring is a resonance ring at the 
UHR or is not a ring at all. 

4.3. Episodic inflov^r to the nuclear disk 

We have shown that the northwest stream is a 
~ 2 X 10^ Mq gas clump on its way from the 3 kpc ring 
to the nuclear disk. At an infall speed of 200 km s~^ 
and a distance of 1 kpc from the center, most of this 
mass should soon collide and merge with the nuclear 
disk; part of the stream should already be colliding 
with the outskirts of the disk. A crude estimate of 
the rate of accretion to the disk during this collision 
is M = M/{l/V) ^ 50 Mq yr"\ where I is the size 
of the northwest stream and we assume that most 
of the stream material will merge with the nuclear 
disk during its first passage of pericenter ( |4.l| ). The 
absence of inflowing gas along the full length of the 
bar means that the high accretion rate lasts for only 
a few Myr, and that inflow must be largely episodic 



rather than continuous. The next episode after the 
present one may well occur as gas in the possible "east 
stream" moves inward. While a steady supply of ten- 
uous gas may have escaped detection, the impending 
increase of the nuclear disk's mass by ~ 10 % seems 
likely to be more significant. As this fueling event 
progresses, we would expect it to produce cloud col- 
lisions, shocks, and gravitational instabilities in the 
nuclear disk, enhancing the star formation rate and 
(perhaps) the rate of inflow to the active nucleus. If 
the supply of gas to the nucleus is similarly episodic 
in other galaxies, we have a natural mechanism to 
regulate the intermittent triggering of nuclear star- 
bursts, which Balzano (1983) infers have a duty cycle 
of ~ 5 %. A connection between kpc-scale fueling 
and nuclear activity is likely to be more indirect at 
best, since accretion luminosity can be significantly 
boosted by the disruption of an individual star (e.g. 
Eracleous et al. 1995) or gas cloud (e.g. Bottema & 
Sanders 1986). 

5. Conclusions 

We have observed CO (1-0) emission from the LINER 
galaxy NGC 5005 to study the gas dynamics in the 
central 6kpc (1'). 

1. In its central 6 kpc, NGC 5005 has three distinct 
molecular gas components — a roughly circular 
ring of radius ^ 3 kpc and gas mass '--^ 2 x 10^ 
Mq; a central condensation with radius ~ 1 kpc 
and gas mass ~ 2 x 10^ Mq; and a component 
with size ~ 0.5 kpc and gas mass ~ 2 x 10^ Mq, 
located about 6" northwest of the nucleus. The 
nuclear disk shows weaker CO(l-O) emission in 
its innermost 100 - 200 pc. 

2. The 3 kpc ring shows noncircular velocities, with 
50 km s^^ inward motions observed on the mi- 
nor axis. The northwest stream has a very large 
noncircular motion; its line-of-sight velocity is 
~ 150 km s~^ different from that of the circu- 
lar galactic rotation. The nuclear disk has a 
linewidth ~ 30 % larger than that of the 3 kpc 
CO ring and of the entire galaxy's HI emission. 
Noncircular motions are also evident in the first 
moment map of the disk. 

3. A stellar bar of length ~ 5 kpc may account 
for these gas distributions and their kinematics. 
The bar would run roughly east- west, placing 
the northwest stream and an associated straight 



dust lane at its leading edge. The nuclear disk is 
mainly formed from X2 orbits elongated perpen- 
dicular to the bar. The disk's large linewidth is 
probably due to clouds on elliptical orbits which 
move along our line of sight as they reach peri- 
center, rather than (necessarily) to a large mass 
concentration at the galactic center. The north- 
west stream lies on an xi orbit elongated par- 
allel to the bar; it is falling towards the nuclear 
disk, with which it will ultimately collide and 
merge. 

4. The 3 kpc ring may form at the ultraharmonic 
resonance of the bar or (rather improbably) con- 
stitute a pair of spiral arms which originate 
at its ends. Noncircular motion would be ex- 
plained in these scenarios by orbital motions in 
an elongated ring and by spiral-arm streaming, 
respectively. 

5. The supply of gas to the nuclear disk is episodic 
rather than continuous. Gas from the outer 
galaxy will accumulate in the 3 kpc ring and oc- 
casionally fall to the nuclear disk, as in the case 
of the northwest stream, due to the action of the 
bar. The present fueling event will ultimately 
increase the mass of the nuclear disk (by about 
10 %) and impact star formation there within 
several Myr. A subsequent fueling event may 
occur if gas in the possible "east stream" is truly 
flowing along the bar. Such episodic supply of 
gas may be responsible for triggering intermit- 
tent nuclear starbursts in other disk galaxies. 

A future paper (Baker et al. 1999) will present fur- 
ther analysis of the nuclear disk and its relation to 
fueling of the LINER, based on 1" resolution maps of 
multiple molecular lines. 
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Fig. 1. — CO (1-0) zeroth and first moment maps of 
NGC 5005. Syntfiesized beams are sliown in tlie lower 
left corners; large crosses mark the dynamical center 
at ai95o = 13'^08'"37;67 and 61950 = +37°19'28'.'34. 
(Top row) Low-resolution (4") maps. The ellipse is 
centered on the dynamical center, and has major and 
minor axes of 60" and 60" cos 62° with a position an- 
gle of 65°, which corresponds to a circle in the galac- 
tic plane after deprojection. (Middle row) Medium- 
resolution (2"7) maps. (Bottom row) High-resolution 
(1'.'5) maps of the nuclear gas disk. The small crosses 
mark the mean position of the radio continuum nu- 
cleus. 

Fig. 2. — Low-resolution channel maps of CO(l-O) 
emission from NGC 5005. Each panel is labeled with 
the LSR velocity in km s^"'^. The synthesized beam is 
shown in the lower left corner of the first panel. The 
ellipse is identical to that shown in Figure nl 

Fig. 3.— Optical image of NGC 5005 (A = 550 nm). 
North is to the top and east is to the left. 

Fig. 4. — The low-resolution CO(l-O) map overlaid 
on the HST ^-band image of NGC 5005. Dark dust 
lanes and molecular gas correlate well. 

Fig. 5. — Position-velocity cuts along the major and 
minor axes of the galaxy (p. a. 65° and 155°, respec- 
tively), made from the medium-resolution data. Po- 
sition is measured from the dynamical center, and 
velocity is relative to 932.8 km s~^. 

Fig. 6.— The central 14" x 14" of NGC 5005. The 
high-resolution CO(l-O) map is overlaid on the HST 
y-band image. The cross marks the radio nucleus. 

Fig. 7. — Medium-resolution CO (1-0) channel maps 
for channels containing emission from the northwest 
stream, the centroid of which is shown as a cross. 
Each panel is labeled with the LSR velocity in km s~^. 
The synthesized beam is shown in the lower left cor- 
ner of the first panel. The ellipse is identical to that 
shown in Figure n^. 

Fig. 8. — Schematic view of the gas distributions 
in NGC 5005. The 3kpc ring, the nuclear molecular 
disk, and the northwest stream are seen in CO, as is 
a possible "east stream." A dust lane between the 
nuclear disk and the 3kpc ring (coincident with the 
northwest stream) and the outer spiral arms are sug- 
gested from dark lanes in optical images. A bar or 



oval distortion shown as a dashed ellipse is inferred 
from the gas distribution and kinematics. The major 
axis of the galaxy has p. a. 65°; the northwest side of 
the galaxy is near to us; and the sense of rotation is 
counterclockwise. 

Fig. 9. — Position-velocity cut along the bar (p. a. 
80°) through the dynamical center, made from the 
low-resolution data. The cut is 5" wide and contains 
the northwest stream at about (-1-5", +100 km s~^) 
and the possible "east stream" at about (—20", -1-100 
km s^^). The velocity zero-point is 937.9 km s^^. 

Fig. 10. — (Left) Model orbits in a barred poten- 
tial, in a frame rotating with the bar. The bar 
runs horizontally (along the dashed line) and ro- 
tates counterclockwise. The potential is $(r, 9) = 
— (1 + £cos20)log(l + r^)/2, where {r,9) are polar 
coordinates rotating with the bar and the strength of 
the bar is £ = 0.04. The pattern speed of the bar is 
0.1, which places the inner and outer ILR, UHR, and 
corotation radii at 1.1, 2.1, 6.3, and 10, respectively. 
The orbits are calculated using a damped orbit model 
(Wada 1994). Horizontally elongated orbits are xi 
orbits, while vertically elongated ones are X2 orbits. 
(Right) Position-velocity cut (of full spatial width 8) 
through the model galaxy at left, along the bar major 
axis. The solid line shows the circular velocity curve. 
Gas particles on different orbits are plotted with dif- 
ferent symbols: crosses are on X2 orbits, diamonds are 
on xi orbits at the leading edge of the bar, and dots 
are on xi orbits at the traihng edge of the bar. Gas 
on X2 orbits and on xi orbits at the leading edge of 
the bar (i.e. crosses and diamonds) form a "tilted X" 
pattern like that seen in Figure || for NGC 5005. Ve- 
locities have been corrected for inclination, assuming 
that the line of nodes and the bar major axis differ 
by 30° in the plane of the galaxy. 

Fig. 11. — Face-on view of NGC 5005, made from the 
low-resolution map. Deprojection assumed i — 62° 
and P. A. = 65°. Primary-beam correction has been 
applied. In this figure, the line of nodes is horizontal, 
the reference circle has a 30" radius, and the line at 
p. a. « —45° shows the approximate position angle of 
the bar (east- west in the sky plane) . The approaching 
side of the galaxy is to the right, and the near side is 
at top. Gas rotates counterclockwise. 

Fig. 12. — Velocities on the 3kpc ring are plotted as 
a function of azimuthal angle in the galactic plane. 
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where 6 is measured counterclockwise from the re- 
ceding major axis (p. a. — 65°). Data points (as- 
terisks) are mean hne-of-sight velocities for the ring 
{R = 25" - 35") with ±1(T standard deviations, mea- 
sured from the low-resolution CO data. The solid 
curve shows velocities expected for circular rotation 
with azimuthal velocity V^z = 295 km s^^ in the 
galactic plane (T^ys = 952 km s~^). The dotted curve 
shows the expectation for an oval orbit with Vrad = 
-60sin2(6l - 30°) and Kz = 295 - 30cos2(6' - 30°). 
At bottom are the residual velocities for each model. 
The 30° phase offset corresponds to the difference be- 
tween line of nodes and bar major axis which was used 
to correct Figure Eo^ for inclination. 
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Table 1 






Parameters of NGC 5005 




Parameter 


Value 


Ref. 


Hubble type 


SAB(rs)bc 


(1) 


Distance [Mpc] 


21.3 


(2) 


P.A. (isophotal) [°] 


65 


(1) 


inclination (isophotal) [°] 


62 


(1) 


Vsys [km s-1]- 


952 


(1) 


5° [mag] 


10.19 


(1) 


^CO(i-o)(45") [Jykms-i] 


697 ± 90 


(3) 


^co(i-o)(23") [Jykms-i] 


435 ± 31 


(4) 



Note. — References are (1) dc Vaucouleurs 
et al. (1991); (2) TuUy et al. (1988); (3) Young 
et al. (1995), 45" (FHWM) beam centered at 
ai95o = 13'^08'"37!6 and (5i95o = +37°19'25"; (4) 
Brainc et al. (1993), 23" (FWHM) beam centered 
at ai95o = 13h08'"37!7 and ^1950 = +37°19'28"5. 

''HI velocity, converted to the LSR velocity in 
the radio convention. 



Table 2 
Parameters of CO(l-O) maps 



Resolution 


Configs 


& Weighting 




Beam 




Kes 




a 




/45" 


/23" 


(1) 




(2) 






(3) 




(4) 




(5) 




(6) 


(7) 










II 


II 





km s-i 


mjy bm 


-1 


niK 


% 


% 


Low 


L+H; 


NA 


4.57 


3.46 


-21.25 


31.1 


9.9 




58 


75 


92 


Medium 


L+H+U; 


UN. 


robust 2.0 


2.88 


2.44 


-18.88 


10.4 


13.8 




181 


73 


94 


High 


L+H+U; 


UN, 


robust -0.5 


1.55 


1.42 


-89.95 


20.7 


16.8 




705 


47 


75 



Note. — Col. (2): L,H,U are the observing configurations. UN = uniform weighting, NA = natural weighting. Col. 
(3): Full width at half maximum along the major and minor axes of the beam, and position angle of the major axis. Col. 
(4): Velocity width of each channel map. Col. (5): \a noise level expressed in intensity and brightness temperature. 
Col. (6): Fraction of flux recovered in the data cube, out of total seen by the FCRAO 14 m telescope in a 45" FWHM 
beam (Young et al. 1995). Col. (7): Fraction of flux recovered in the data cube, out of total seen by the IRAM 30 m 
telescope in a 23" FWHM beam (Braine et al. 1993). 
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Table 3 
Position of the nucleus 



# 


R.A. (B1950) 


Dec, 


. (B1950) 


Resolution 


Aobs 


Ref. 




h 


in 


S 


O 


/ 


II 


// 


cm 




1 


13 


08 


37.68 


+37 


19 


2%.ll 


1 


6 


(a) 


2 


13 


08 


37.67 


+37 


19 


28.4 


1.5 


20 


(b) 


3 


13 


08 


37.69 


+37 


19 


28.6 


1.0 


6,20 


(c) 


4 


13 


08 


37.651 


+37 


19 


28.80 


5 


20 


(d) 


5 


13 


08 


37.673 


+37 


19 


28.64 






(e) 


6 


13 


08 


37.67 


+37 


19 


28.34 


1.5 


0.3 


(f) 



Note. — Sources are (a) van der Hulst et al. (1981); (b) Condon et al. 
(1990); (c) mean of 6 and 20 cm positions in Vila et al. (1990); (d) 
VLA/FIRST survey (White et al. 1997), converted from J2000 to B1950; 
(e) arithmetic mean of 1,2,3, and 4; (f) this paper: dynamical center de- 
termined from the high-resolution CO (1-0) data. 



Table 4 
Kinematic parameters of the nuclear disk 



Parameter 


Value 


Dynamical center aigso 


13'^08™37?67 


(5l950 


+37°19'28'.'34 


P.A. [°] 


70 


inclination [°] 


53 


V^ys [km s-1] 


928^ 


Velocity range [km s~^] 


565 - 1301*^ 


Ay [km s-i] 


736^ 



Note. — All values except the velocity 
range and the linewidth are determined 
from a fit to the high-resolution first mo- 
ment map in Figure nl The northwest 
stream is excluded from the fit. For the re- 
mainder of our analysis, we adopt the very 
similar isophotal position angle and incli- 
nation in Table H. 

'^LSR velocity defined in the radio con- 
vention. 

''Includes channels with emission at the 
3(7 level in the medium-resolution data 
cube. 
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Table 5 
Gas and dynamical masses 


Radius'* M^yn^ 
" kpc IQio Mq 


10^ M© 


Mgas/Afdyn 
% 


5 0.5 2.1 
30 3.1 6.3 


1.4 
4.0 


7 
6 



'^Radius in the galactic plane. 

Rotation velocities are measured from the ma- 
jor axis position-velocity diagram in Figure ^ The 
dynamical mass in the central 5" is probably over- 
estimated because of noncircular motion (see jj4.l|). 



'^CO fluxes are measured from the medium- 
resolution data for R < 5" and from the low- 
resolution data for R < 30". 



16 



I I I I 








"P_ 



J I I I I I I 



I I I 




"^ 



J I I I I I L 



^^^Q ^^^^ ^Q^^ ^^^Q 



■ lib, llb,IBb, I^U 




IIIUI 31 K 3- 

rHHTH^aiBDIIHI 



■ M.^Lfl.JULJHh 



1 1 1 1 ^1 

-| 1 1 1 1 1 1 1 



O 




J I I I I I L 



- i 



8 



I I !:!!!?Sg 
T 1 1 1 1 1 1 1 




■_b X-. 1. ^^^ 



x-j a.-, a III 



iM. iM i-j i-.a r^ 



Lhi ■■ LUBIBiD 



17 

Fig. 1. — Quality of this figure was lowered for astro-ph. A PS preprint with the full-resolution figures is avaiable 
from http://www.ovro.caltech.edu/mm/science/science.html . 
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18 
Fig. 2. — Quality of tliis figure was lowered for astro-ph. A PS preprint with the full-resolution figures is avaiable 
from http: //www. ovro . caltech.edu/mm/science/science .html . 




Fig. 3. — Quality of this figure was lowered for astro-ph. A PS pre print with the full-resolution figures is avaiable 



from http: //www. ovro . caltech.edu/mm/science/science .html 
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Fig. 4. — Quality of this figure was lowered for astro-ph. A PS pre print with the full-resolution figures is avaiable 
from http: //www. ovro . caltech.edu/mm/science/science .html . 
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Fig. 5. — Quality of this figure was lowered for astro-ph. A PS preprint with the full-resolution figures is avaiable 
from http://www.ovro.caltech.edu/min/science/science.html . 
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Fig. 7. — Quality of this figure was lowered for astro-ph. A PS pre print with the full-resolution figures is avaiable 
from http : //www . ovro . caltech . edu/mm/science/science . html . 
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Fig. 11. — Quality of this figure was lowered for astro-ph. A PS preprint with the full-resolution figures is avaiable 
from http://www.ovro.caltech.edu/inin/science/science.html . 
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